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Abstract. 

The first simultaneous operation of the AURIGA detector and the LIGO 
observatory was an opportunity to explore real data, joint analysis methods 
between two very different types of gravitational wave detectors: resonant bars 
and interferometers. This paper describes a coincident gravitational wave burst 
search, where data from the LIGO interferometers are cross-correlated at the time 
of AURIGA candidate events to identify coherent transients. The analysis pipeline 
is tuned with two thresholds, on the signal-to-noise ratio of AURIGA candidate 
events and on the significance of the cross-correlation test in LIGO. The false 
alarm rate is estimated by introducing time shifts between data sets and the 
network detection efficiency is measured with simulated signals with power in the 
narrower AURIGA band. In the absence of a detection, we discuss how to set 
an upper limit on the rate of gravitational waves and to interpret it according to 
different source models. Due to the short amount of analyzed data and to the high 
rate of non-Gaussian transients in the detectors noise at the time, the relevance 
of this study is methodological: this was the first joint search for gravitational 
wave bursts among detectors with such different spectral sensitivity and the first 
opportunity for the resonant and interferomctric communities to unify languages 
and techniques in the pursuit of their common goal. 



1. Introduction 

Gravitational wave bursts are short duration perturbations of the space-time metric 
due to such catastrophic astrophysical events as supernova core collapses [1] or the 
merger and ringdown phases of binary black hole coalescences [2, 3]. Over the past 
decade, the search for these signals has been independently performed by individual 
detectors or by homogeneous networks of resonant bars [4] or laser interferometers [5, 
6, 7, 8, 9]. The first coincident burst analysis between interferometers with different 
broadband sensitivity and orientation was performed by the TAMA and LIGO 
Scientific Collaborations [10]. That analysis required coincident detection of power 
excesses in at least two LIGO interferometers and in the TAMA detector in the 700- 
2000 Hz frequency band, where all sensitivities were comparable. The upper limit 
result accounted for the different antenna patterns with a Monte Carlo estimate of 
detection efficiency for sources uniformly distributed in the sky. 

This paper describes a joint burst search in a more heterogeneous network, 
comprised of LIGO and AURIGA. Although this is the first joint search with a 
resonant antenna, bar data has been cross-correlated with LIGO data in the search 
for gravitational wave stochastic background, with the ALLEGRO detector. That 
search provided what's, to date, the most competitive stochastic upper limit in the 
905-925 Hz frequency band [14]. 

LIGO consists of three interferometers, two co-located in Hanford, WA, with 2 km 
and 4 km baselines and one in Livingston, LA, with a 4 km baseline, sensitive between 
60 and 4000 Hz with best performance in a 100 Hz band around 150 Hz. AURIGA 
is a bar detector equipped with a capacitive resonant transducer, located in Legnaro 
(PD), Italy. In 2003 the AURIGA detector resumed data acquisition after upgrades 
that enlarged its sensitive band to 850-950 Hz, from the ~2 Hz bandwidth of the 
1997-1999 run [11, 12, 13]. 

Due to the different spectral shapes, an interferometer-bar coincident search is 
only sensitive to signals with power in the bar's narrower band. The LIGO-AURIGA 
analysis thus focused on short duration (<20 ms) transients in the 850-950 Hz band, 
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with potential target sources like black hole ringdowns [2] and binary black hole 
mergers [15, 16]. 

Another important difference between bars and interferometers is the sky 
coverage, which depends on the detectors' shape and orientation. Figure 1 shows the 

antenna pattern magnitude + of the AURIGA and LIGO-Hanford (LHO) 

detectors, as a function of latitude and longitude. Since directions of maximum LIGO 
sensitivity overlap with the larger portion of the sky visible to AURIGA, a coincident 
search is not penalized by differences in antenna pattern. However, adding AURIGA 
to the detector network does not improve its overall sky coverage either, due to the 
AURIGA sensitivity, which is ~ 3 times worse than LIGO [17, 18]. 

Despite the different sensitivity and bandwidth, a coincident analysis between 
LIGO and AURIGA has the potential to suppress false alarms in the LIGO network, 
thus increasing the confidence in the detection of loud signals and making source 
localization possible, with triangulation. Collaborative searches also increase the 
amount of observation time with three or more operating detectors. For this reason, 
and to bring together the expertise of two traditions in burst analysis, the AURIGA 
and LIGO Scientific Collaborations pursued a joint search. 

The analysis described in this paper follows the all-sky approach described in [18], 
where data from two or three LIGO interferometers are cross-correlated at the time of 
AURIGA candidate events. This method was tested on data from the first AURIGA 
and LIGO coincident run, a 389 hour period between December 24, 2003 and January 
9, 2004, during the third LIGO science run S3 [7] and the first run of the upgraded 
AURIGA detector [12, 13]. Only a portion of this data was used in the joint burst 
search, because of the detectors' duty factors and the selection of validated data 
segments which was independently performed by the two collaborations [17, 19]. The 
effective livetime available for the analysis was: 

- 36 hours of 4-fold coincidence between AURIGA and the three LIGO 
interferometers; 

- 74 hours of 3-fold coincidence between AURIGA and the two LIGO Hanford 
interferometers, when data from the LIGO Livingston detector was not available. 

Other three-detector combinations including AURIGA were not considered, due to the 
low duty factor of the LIGO Livingston interferometer in S3. The 4-fold and 3-fold 
data sets were separately analyzed and the outcome was combined into a single result. 

Figure 2 shows the best single-sided sensitivity spectra for LIGO and AURIGA 
in the 800-1000 Hz band at the time of the coincident run. The AURIGA spectrum 
contained spurious lines, due to the up-conversion of low frequency seismic noise. 
These lines were non stationary and could not always be filtered by the AURIGA 
data analysis; for this reason, a large portion of the data (up to 42%) had to be 
excluded from the analysis, with significant impact on the livetime [17, 19]. The 
largest peak in each LIGO spectrum is a calibration line, filtered in the analysis. The 
amount of available LIGO livetime was limited by several data quality factors, such 
as data acquisition problems, excessive dust at the optical tables, and fluctuations of 
the light stored in the cavities, as described in [7]. 

Due to the short duration of the coincidence run and the non-optimal detector 
performances, the work described in this paper has a methodological relevance. On 
the other hand, it is worth pointing out that the three-fold coincidence between 
AURIGA and the two Hanford interferometers, when Livingston was offline, allowed 
the exploration of some data that would not have been searched otherwise. 
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2. The analysis pipeline 

The joint analysis followed a statistically blind procedure to avoid biases on the result: 
the pipeline was tested, thresholds were fixed and procedural decisions were made 
before the actual search, according to the following protocol [18, 19]. 

(i) AURIGA provided a list of burst candidates (triggers) in the validated observation 
time. The triggers were identified by matched filtering to a (S-like signal, with 
signal-to-noise ratio threshold SNR> 4.5. Triggers at lower SNR were not 
included in this analysis, since their rate and time uncertainty increased steeply to 
unmanageable levels, with negligible improvement in detection efhciency. Special 
attention was required, in this run, to address non-stationary noise with data 
quality vetoes that were not needed in subsequent runs [17, 19]. The resulting 
events were auto-corrclated up to about 300 s; this effect, particularly evident 
for high SNR events, was due to an imperfect suppression of the non-stationary 
spurious lines on short time scales. 

(ii) Data from the three LIGO interferometers at the time of AURIGA triggers were 
cross-correlated by the r-statistic waveform consistency test [20], a component 
of the LIGO burst analysis [6, 7] performed with the CorrPower code [21]. The 
test compares the broadband linear cross-correlation r between two data streams 
to the normal distribution expected for uncorrelated data and computes its p- 
value, the probability of getting a larger r if no correlation is present, expressed 
as r = —log]^Q(p- value). When more than two streams are involved, F is the 
arithmetic mean of the values for each pair. The cross-correlation was performed 
on 20, 50 and 100 ms integration windows, to allow for different signal durations. 
Since the source direction was unknown, the integration windows were slid around 
each AURIGA trigger by 27 ms+at , sum of the light travel time between AURIGA 
and Hanford and of the estimated la timing error of the AURIGA trigger. The 
value of (Tt depended on the SNR of each trigger, typically in the 5 — 40 ms range, 
with an average value of 17 ms. The resulting F was the maximum amongst all 
time slides and integration windows. Only triggers above the minimal analysis 
threshold of F > 4 were considered as coincidences. 

(iii) A cut was applied on the sign of the correlation between the two Hanford 
interferometers, which must be positive for a gravitational wave signal in the two 
co-located detectors. This cut, also used in the LIGO-only analysis [7], reduced by 
a factor ~2 the number of accidental coincidences, with no effect on the detection 
efRciency. 

(iv) The data analysis pipeline was first applied, for testing purposes, to a playground 
data set [22], which amounted to about 10% of the livetimc and was later excluded 
from the data set used in the analysis. 

(v) The false alarm statistics were estimated on off-source data sets obtained by time 
shifting the LIGO data; more details are provided in section 2.1. 

(vi) The detection efficiency was estimated with Monte Carlo methods [19], by 
simulating a population of sources uniformly distributed in the sky and in 
polarization angle, as described in section 2.2. 

(vii) The analysis timing consisted of setting two thresholds: on the SNR of the 
AURIGA candidate events and on the LIGO F value. Details on the tuning 
procedure are available in section 2.3. 
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(viii) The statistical analysis plan was defined a priori, with decisions on which 
combination of detectors to analyze (4-fold and 3-fold) and how to merge the 
results, the confidence level for the null hypothesis test, and the procedure to 
build the confidence belt, as described in section 2.5. 

(ix) Once analysis procedure and thresholds were fixed, the search for gravitational 
wave bursts was applied to the on-source data set. The statistical analysis led to 
confidence intervals which were interpreted as rate upper limit versus amplitude 
curves. A posteriori investigations were performed on the on-source results (see 
section 3.1), but these follow-up studies did not affect the statistical significance 
of the a priori analysis. 

2.1. Accidental coincidences 

The statistics of accidental coincidences were studied on independent off-source data 
sets, obtained with unphysical time shifts between data from the Livingston and 
Hanford LIGO detectors and AURIGA. The two Hanford detectors were not shifted 
relative to each other, to account for local Hanford correlated noise. The shifts applied 
to each LIGO site were randomly chosen between 7 and 100 seconds, with a minimum 
separation of 1 second between shifts. Hanford-Livingston shifts in the 4-detector 
search also had to differ by more than 1 second. The livetime in each shifted set 
varied by a few percent due to the changing combination of data quality cuts in the 
various detectors. The net live time used in the accidental rate estimate was 2476.4 h 
from 74 shifts in the four-detector search and 4752.3 h from 67 shifts in the three- 
detector search. 

Figure 3 shows scatter plots of the LIGO F versus the AURIGA SNR for 
background events surviving the cut on the Hanford-Hanford correlation sign, in the 
4-detector and in the 3-detector configurations. The regions at SNR < 4.5 and F < 4 
are shaded, as they are below the minimal analysis threshold. 

The number of off-source accidental coincidences in each time shift should 
be Poisson distributed if the time slide measurements are independent from each 
other. For quadruple and for triple coincidences, a test compared the measured 
distributions of the number of accidentals to the Poisson model. The test included 
accidental coincidences with F > 4 and F > 7.5 for 4-fold and 3-fold coincidences, 
respectively. These thresholds were lower than what was used in the coincidence 
search (sec. 2.3), to ensure a sufficiently large data sample, while the AURIGA 
threshold remained at SNR > 4.5. The corresponding p-values were 34% and 6.5%, 
not inconsistent with the Poisson model for the expected number of accidentals. 

2.2. Network detection efficiency 

The detection efficiency was estimated by adding software-generated signals to real 
data, according to the LIGO Mock Data Challenge procedure [23]. The simulation 
generated gravitational waves from sources isotropically distributed in the sky, with 
azimuthal coordinate uniform in [0, 27r], cosine of the polar sky coordinate uniform 
in [—1, 1] and wave polarization angle uniform in [0, tt]. Three waveform classes were 
considered [17, 19]: 

(i) Gaussians with linear polarization: 

I /i+(t) = Va;=e-(*-*«)V-^ 
I h^{t) = Q 
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with T — 0.2 ms. 
(ii) sine-Gaussians with linear polarization: 

h+it) = hpeak e-(*-*«)V-' sin(2^/o(t - to)) 



h^{t) = 

with fo = 900 Hz, r = 2//o = 2.2 ms and Q = v^tt/ot = 8.9. In this analysis we 
also tested cosine-Gaussians (with the sin replaced by cos), and found the same 
sensitivity as for sine-Gaussian waveforms. 

(iii) Damped sinusoids with circular polarization: 

^ f Vafciif^ cos[2^/o(i-to)]e-(*-*°)/^ i-to>0, 

h+[t) = < 1 2 

[ Vafe i + f^ C0s[2^/o(t - to)] clO(*-*o)/- t - to < , 

, (..^{ hp,ak COS i sin[27r/o(t - to)] e-(*-*«)/" t - to > 
^ ^ ^ \ /ipeafc COS L sin[27r/o(t - to)] eio(*-*»)/- t - to < 

with fo = 930 Hz, r ~ 6 ms and cost uniformly distributed in [—1, 1], t being the 
inclination of the source with respect to the line of sight. 

Although no known astrophysical source is associated with Gaussian and sine-Gaussian 
waveforms, they are useful because of their simple spectral interpretation and they 
are standard test waveforms in LIGO burst searches. Damped sinusoids are closer to 
physical templates [2, 15, 16]. 

The signal generation was performed with the LIGO LDAS software [24]; the 
waveforms were added to calibrated LIGO and AURIGA data and the result was 
analyzed by the same pipeline used in the search. For each waveform class, the 
simulation was repeated at different signal amplitudes to measure the efficiency of 
the network as a function of the square root of the burst energy: 



Kss = J dt [hlit) + hi it)] ^ J 2 / df h\{f) + hi (/) 




2.3. Analysis tuning 

The analysis thresholds were chosen to maximize the detection efficiency with an 
expected number of accidental coincidences smaller than 0.1 in each of the three and 
four detector searches. Figure 4 shows contour plots of the number of accidental 
coincidences expected in the on-source data set, the original un-shifted data that may 
include a gravitational wave signal, as a function of the F and SNR thresholds. This 
quantity is the number of accidental coincidences found in the time shifted data, scaled 
by the ratio of on-source to off-source livetimes. The plots also show the detectability 
of sine-Gaussian waveforms, expressed as /irss50%, the signal amplitude with 50% 
detection probability. 

For all tested waveforms, the detection efficiency in the 4-fold and 3-fold searches 
are the same, within 10%; their value is dominated by the AURIGA sensitivity at 
SNR > 4.5. This observation, together with the shape of the accidental rate contour 
plots, indicates that the best strategy for the suppression of accidental coincidences 
with minimal impact on detection efficiency is to increase the F threshold and leave 
the SNR threshold at the exchange value of 4.5. The analysis thresholds were chosen 
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to yield the same accidental rate in the two data sets: F > 6 and SNR > 4.5 for the 
4-fold search and T > 9 and SNR > 4.5 for the 3-fold search. 

It was decided a priori to quote a single result for the entire observation time 
by merging the 4-fold and the 3-fold periods. The number of expected accidental 
coincidences in the combined on-source data set is 0.24 events in 110.0 hours, with a 
Icr statistical uncertainty of 0.06. Detection efficiencies with the chosen thresholds are 
listed in table 1. 



Waveform 




[10-20HZ-1/2] 


hrss90% [10 


~19Hz-1/2] 




4-fold 


3-fold 


4-fold 


3-fold 


sine-Gaussians 


5.6 


5.8 


4.9 


5.3 


Gaussians 


15 


15 


10 


11 


damped sinusoids 


5.7 


5.7 


3.3 


3.4 



Table 1. Signal amplitudes with 50% and 90% detection efficiency in the 4-fold 
and 3-fold searches at the chosen analysis thresholds of SNR > 4.5 and F > 6 
(4- fold) or r > 9 (3- fold). These numbers do not include a ~ 10% systematic 
error, due to calibration uncertainties, and a ~ 3% Icr statistical error, due to the 
small number of simulated signals. 



2.4- Error propagation 

The detection efficiency, in a coincidence analysis, is dominated by the least sensitive 
detector, in this case AURIGA. Monte Carlo efficiency studies show that only a small 
fraction of the simulated events are lost due to the LIGO F threshold; these events 
were missed because their sky location and polarization were in an unfavorable part of 
LIGO's antenna pattern. Since most of the simulated events were cut by the AURIGA 
threshold, the main source of systematics in this analysis is the calibration error in 
AURIGA data, estimated to be 10%. 

In addition, there is a statistical error due to the simulation statistics and to the 
uncertainty on the asymptotic number of injections after the veto implementation. 
The 1(7 statistical error on the numbers in table 1 is about 3%. 

Both systematic and statistical errors were taken into account in the final 
exclusion curve in figure 7. The systematic error is propagated from the fit of the 
efficiency curve to a 4-parameter sigmoid [6, 7, 8]. The fit parameters were worsened 
to ensure a 90% confidence level in the fit, following the prescriptions in [27]. An 
additional, conservative shift to the left was applied to account for the 10% error on 
the calibration uncertainty, which is the dominant error. 

2.5. Statistical interpretation plan 

In compliance with the blind analysis approach, the statistical interpretation was 
established a priori. The procedure is based on a null hypothesis test to verify that 
the number of on-source coincidences is consistent with the expected distribution of 
accidentals, a Poisson with mean 0.24. We require a 99% test significance, which 
implies the null hypothesis is rejected if at least 3 coincidences are found. 
The set of alternative hypotheses is modeled by a Poisson distribution: 

P(n|A0 = (/i + &)"exp[-(A* + 5)]/,i! (1) 
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where the unknown ^ is the mean number of counts in excess of the accidental 
coincidences, which could be due to gravitational waves, to environmental couplings 
or to instrumental artifacts. Confidence intervals are established by the Feldman and 
Cousins method with 90% coverage [25]. Uncertainties on the estimated accidental 
coincidence number are accounted for by taking the union of the two confidence belts 
with b = 0.24 ± 3cr, with a = 0.06 . 

The confidence belt was modified to control the false alarm probability according 
to the prescription of the null hypothesis test: if less than 3 events are found, and the 
null hypothesis is confirmed at 99% C.L., we accept the upper bound of the Feldman 
and Cousins construction but we extend its lower bound to regardless of the belt 
value. The resulting confidence belt, shown in figure 5, is slightly more conservative 
than the standard Feldman and Cousins belt for small values of the signal /i. The 
advantage of this modification is to separate the questions of what is an acceptable false 
detection probability and what is the required minimum coverage of the confidence 
intervals [26]. 

An excess of on-source coincidences could be due to various sources, including 
instrumental and environmental correlations; the rejection of the null hypothesis or a 
confidence interval on fi detaching from zero do not automatically imply a gravitational 
wave detection. A detection claim requires careful follow-up studies, to rule out all 
known sources of foreground, or independent evidence to support the astrophysical 
origin of the signal. On the other hand, an upper limit on fj, can be interpreted as 
an upper limit on the number of GWs; therefore the upper bound of the confidence 
interval can be used to construct exclusion curves. 



3. Results 



The final step consists of analyzing the on-source data sets. No gravitational-wave 
candidates were found in this search, consistent with the null hypothesis. The resulting 
90% CL upper limit is 2.4 events in the on-source data set. or 0.52 events/day in the 
combined 3-fold and 4-fold data sets. 

Figure 6 shows the combined efficiency for this search as a function of the signal 
amplitude for the waveforms described in section 2.2, a weighted average of the 
detection efficiency of 3-fold and 4- fold searches: 

/, N e(hrss)4fold • T4fold + £(hi.ss)3fold • Tafold 

£(,nrssj4fold+3fold — 7^ (^j 

-L4fold + -Lsfold 

The 90% C.L. rate upper limit, divided by the amplitude-dependent efficiency, 
yields upper limit exclusion curves similar to those obtained in previous searches [6, 10]. 
Figure 7 compares the sine-Gaussian exclusion curves found in this search to those 
from S2 in LIGO and LIGO-TAMA. The waveform used here peaks at 900 Hz, while 
the previous searches used a sine-Gaussian with 850 Hz central frequency. We verified 
analytically that the AURIGA detection efficiencies for Q=9 sine-Gaussians at 850 Hz 
and 900 Hz agree within 10%; no significant difference is to be expected for the large 
band detectors. 

The asymptotic upper limit for large amplitude signals is inversely proportional 
to the observation time. The value for this search with 90% C.L. is 0.52 events/day, 
to be compared to 0.26 events/day in the LIGO S2 search [6] and 0.12 events/day 
in the LIGO-TAMA search [10]. The lowest asymptotic value was previously set by 
IGEC: 4 X 10""^ events/day, thanks to their longer observation time [4]. 
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The detection efficiency in this search is comparable to the LIGO-only S2 one, 
and a factor 2 worse than the LIGO-only S3 sensitivity. In the lower amplitude region, 
this search is an improvement over the IGEC search, since the AURIGA amplitude 
sensitivity during LIGO S3 was about 3 times better than the typical bar sensitivity in 
the IGEC 1997-2000 campaign (a direct comparison is not possible since IGEC results 
are not interpreted in terms of a source population model). More recent data yielded 
significant improvements in sensitivities, by a factor ~ 10 for the LIGO S4 run [8] and 
a factor - 3 for IGEC-2 [28]. 

3.1. Diagnostics of on-source and off-source data sets 

The agreement between on-source and off-source coincidences was tested comparing 
the r distributions in Figure 8 above the minimal exchange threshold F > 4 and below 
the network analysis threshold (F > 6 for 4-fold and F > 9 for 3-fold) . This a posteriori 
test did not find a disagreement between on-source and off-source distributions. There 
were no 4-fold, on-source events with F > 4. For 3-fold events, the agreement between 
zero-lag and accidental distributions can be confirmed with a Kolmogorov-Smirnov 
test that uses the empirical distribution of accidentals as a model, giving a 0.6 p- 
value. 

In addition, we addressed the question of whether on-source events (foreground) 
modified the distribution of accidentals (background) and biased our estimate. This 
is an issue in the 3-fold AU-H1-H2 analysis where only HI and 112 are cross-correlated 
and the measured background distribution includes instances of Hanford foreground 
events in accidental coincidence with an AURIGA shifted event. As a result, the time 
shift method overestimates the number of accidentals. In this search, however, this 
systematic effect turned out being negligible: the removal of all background events in 
accidental coincidence with on-source 3-fold events with F > 5.5 and SNR> 4.5, did 
not significantly affect the F histogram above threshold. 

The same question could be posed in a different way: how would the background 
histogram change if we had an actual gravitational wave event, with large F? On 
average, the same II1-II2 event appears in ~ 9 background coincidences. A loud 
gravitational wave event, with F above the noise, say F = 12, would have appeared in 
the background histogram 15-20 times as a F peak with a tail at the same Hanford 
time. Such an event would not have been missed, but it would have been noticed in 
the tuning stages. The most significant consequence is that the 3-fold search is not 
truly blind, since a loud signal would easily manifest itself in the tuning data set. 

4. Conclusions 

This paper describes the first joint search for gravitational wave bursts with a 
hybrid network composed of a narrow band resonant bar detector and broadband 
interferometers. This was a rare and valuable opportunity to bring together the 
expertise of the AURIGA and LSC collaborations and explore common methods on 
real data. The addition of the AURIGA detector to the LIGO observatory allowed 
to extend the time coverage of the observations by including also the time periods 
when only two of the three LIGO detectors were operating simultaneously with 
AURIGA (AU-H1-H2). This was possible thanks to the false alarm rate suppression 
contributed by AURIGA. The detection efficiency of this hybrid network for the tested 
source models was about a factor 2 worse than the LIGO-only efficiency, limited by 
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the AURIGA detector. This cost however, turned out smaller than the amplitude 
sensitivity ratio between AURIGA and LIGO during S3 for the same signal types 
(roughly, a factor 3). 

Due to the short observation time, the relevance of this study is methodological. 
The results have been interpreted in terms of source population models and the final 
upper limits are comparable to those set by LIGO alone in previous observations. 
This joint analysis followed a statistically blind procedure to allow an unbiased 
interpretation of the confidence of the results. In particular, the data analysis plan has 
been fixed a priori and the results are confidence intervals which ensure a minimum 
coverage together with a more stringent requirement on the maximum false detection 
probability. 

Any future joint search for bursts by interferometric and resonant detectors, on 
simultaneous long-term observations, would be informed by the techniques developed 
in the ground-breaking work presented in this paper. 
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Figure 1. Antenna pattern magnitude as a function of the longitude and the 
sine of the latitude. Top: AURIGA; bottom: LIGO-Hanford. 
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Figure 2. Best single-sided sensitivity spectra of AURIGA and the three LIGO 
interferometers during the joint observation. In the AURIGA spectrum, several 
spurious lines arc visible while LIGO spectra present calibration lines at 973 Hz for 
the Hanford detectors (LHO-4 km and LHO-2 km) and 927 Hz for the Livingston 
detector (LLO-4 km). 
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Figure 3. Scatter plots of the LIGO F versus the AURIGA SNR for background 
events surviving the cut on the sign of the H1-H2 correlation, in the 4-detoctor 
(left) and in the 3-detector (right) configurations. The regions at SNR < 4.5 and 
r < 4 are shaded, as they are below the minimal analysis threshold and were not 
used in the tuning. 

The vertical structures at large SNR and F < 4 are an artifact, from the cross- 
correlation of LIGO data at different time shifts around the same loud AURIGA 
event. 




AURIGA SNR threshold AURIGA SNR threshold 

Figure 4. Grayscale plot: efficiency to sine-Gaussians waveforms in terms of 
hrss50%. Left: 4-fold; right: 3-fold observations. 

The contour lines indicate the number of accidental coincidences expected in the 
on-source data set. 
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Figure 5. Confidence belt based on the Feldman and Cousins procedure, with 
90% C.L., unknown Poisson signal mean fi and a Poisson background with 
b = 0.24 ± 3(T. The standard Feldman and Cousins construction (solid lines) is 
modified by fixing the maximum false alarm probability to 1%: the lower bound 
is fixed to for 2 measured events. 




.|Q-21 .|Q-20 .|Q-19 .|Q-18 .|Q-17 

hrss[Hz-^''] 

Figure 6. Average efficiency of detection (4-fold+3-fold) versus hrss for the 
considered waveforms: sine-Gaussians at /o = 900 Hz and r = 2.2 ms, Gaussians 
with r = 0.2 ms, and damped sinusoids with /o = 930 Hz and damping time 
T = 6 ms. All sources have been modelled as uniformly distributed in the sky and 
with random polarizations (see section 2.2). The triangles mark the hrss50% for 
sine-Gaussians achieved in the LIGO-only S2 and S3 searches. 
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Figure 7. Upper limits at 90% C.L. on the gravitational wave rate versus h^-sa 
for sine-Gaussian waveforms in different network analyses: AURIGA-LIGO S3 
(dotted line), LIGO-only S2 (solid line) and TAMA-LfGO S2 (dashed line). The 
sine-Gaussians have central frequency 900 Hz (for AURIGA-LIGO S3) and 850 Hz 
(for LIGO-only S2 and TAMA-LIGO S2). In all cases, Q = 8.9. All sources have 
been modelled as uniformly distributed in the sky and with random polarizations 
(see section 2.2) 
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Figure 8. Comparison of F distributions for four-fold (left) and three-fold 
coincidences (right). The black dots correspond to on-source coincidences. The 
stair-step curve with gray dots is the mean accidental background, estimated from 
off-source coincidences, with its la RMS spread (shaded gray area) and the error 
on the mean (thin error bars). The on-source and off-source distributions are 
in good agreement, within the statistical uncertainty. The solid vertical lines 
correspond to the analysis thresholds (F = 6 and F = 9, respectively). The 
analysis was tuned only on events with F > 4 (dotted vertical lines). 



